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sensitive to the impedance of the environment, which can be tuned in situ. We have observed 
Coulomb blockade of Cooper-pair tunneling and even a region of negative differential resistance, 
when the zero-bias resistance R' of the SQUID arrays is much higher than the quantum resistance 
Rk = h/e 2 « 26 kf2. The negative differential resistance is evidence of the coherent single-Cooper- 
pair tunneling within the theory of current-biased single Josephson junctions. Based on this theory, 
we have calculated the I-V curves numerically in order to compare them with the experimental ones 
at Rq 3> Rk- The numerical calculation agrees with the experiments qualitatively. We also discuss 
the R' dependence of the single- Josephson-j unction I-V curve in terms of the superconductor- 
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PACS numbers: 74.50.+r, 73.23.Hk, 73.40.Gk 
Phys. Rev. B 67, 094505 (2003). 



X 



I. INTRODUCTION 



Small-capacitance superconducting tunnel junctions 
provide an ideal system for studying the interplay be- 
tween quantum mechanically conjugate variables: the 
Josephson-phase difference across the junction and the 
charge on the junction electrode. Their behavior is influ- 
enced by dissipation, such as discrete tunneling of quasi- 
particles, coupling to an electromagnetic environment, 
etc. The simplest and the most fundamental example is 
the single junction. Current-voltage {I-V) characteris- 
tics of single junctions have been the subject of extensive 
theoretical investigations^^ A small capacitance C gives 
a large charging energy e 2 /2C, and thus the theory pre- 
dicts that at low enough temperatures, a superconducting 
small-capacitance tunnel junction behaves like an insu- 
lator so long as the junction is isolated from the elec- 
tromagnetic environment. Experimentally, however, the 
observation of this charging effect (Coulomb blockade) 
in a single junction is not straightforward, because the 
measurement leads attached to the single junction should 
have a high enough impedance to isolate the junction. 1 
For this reason, thin-film resistors were employed for the 
leads, and Coulomb blockade was successfully observed in 
single Josephson junctions. 3 Tunnel-junction arrays were 
also tried for the leads, and an increase of differential 
resistance around V = was reportedA^ 

We use one-dimensional (ID) arrays of dc supercon- 
ducting quantum interference devices (SQUIDs) for the 
leads.' 1 The advantage of this SQUID configuration is 
that the effective impedance of an array can be varied 



in situ by applying an external magnetic field perpen- 
dicular to the SQUID loop. Thus, the zero-bias resis- 
tance of the SQUID arrays at low temperatures can be 
controlled over several orders of magnitude. This phe- 
nomenon has been extensively studied in terms of the 
superconductor-insulator (SI) transitioni&iliiS The sin- 
gle junction in our samples, on the other hand, does not 
have a SQUID configuration, and therefore its param- 
eters are practically independent of the external mag- 
netic field. This enables us to study the same single 
junction in different environments. We show that the 
I-V curve of the single junction is indeed sensitive to 
the state of the environment. Furthermore, we can in- 
duce a transition to a Coulomb blockade of the single 
junction by increasing the zero-bias resistance R' of the 
SQUID arrays. When R' is much higher than the quan- 
tum resistances: Rk = h/e 2 w 26 kf2 for quasiparticles 
and Rq = h/{2e) 2 « 6.5 kf2 for Cooper pairs, single- 
junction I-V curve has a region of negative differential 
resistance. The negative differential resistance is a result 
of coherent single-Cooper-pair tunneling from the view- 
point of the theory for current-biased single Josephson 
junctions whose electromagnetic environment has a suf- 
ficiently high impedancei^ 

In the context of the theory, we calculate the I-V 
curves numerically by following Ref. Illl and compare 
with the experimental ones at R' 3> Rk- The nu- 
merical results are in qualitative agreement with the 
experiments, 12 however, some quantitative discrepancies 
suggest that even at R' » Rk, the single-j unction I-V 
curve would be influenced by the environment. Accord- 
ing to the perturbation theorjii 3 ^^ that deals with the 
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effect of a linear-impedance environment Z(u) on sin- 
gle Josephson junctions, Re[Z(w)], or the dissipation, 
plays a key role. An insulating behavior characterized 
by the Coulomb blockade of Cooper-pair tunneling is 
expected for the single junction when Re[i?(w)] 3> Rq, 
i.e., the dissipation is sufficiently small, because the large 
Re[i?(a;)] suppresses the charge fluctuation. The dissipa- 
tive dynamics of single Josephson junctions has also been 
discussed by studying resistively shunted junctions both 
theoreticallySiiSiii and experimentallyii2*iS As the dissi- 
pation is increased (the shunt resistance is decreased), 
the tunneling of the Josephson phase is suppressed, and 
as a result the Josephson junction exhibits a transition 
from an insulator to a superconductor. The theoreti- 
cal phase diagram^ for this SI transition was supported 
by the experimentsi&i^ qualitatively. The superconduct- 
ing behavior has been observed experimentally in single 
Josephson junctions not only by shunting with a resistor 
but also by biasing with a low- impedance circuit^ i.e., 
by achieving a low-impedance electromagnetic environ- 
ment. Thus, we expect that there should be a SI tran- 
sition driven by the electromagnetic environment. We 
discuss the environment-driven SI transition in our sys- 
tem, where the environment is tunable with the SQUID 
arrays. 

II. THEORY 

The Hamiltonian of a single Josephson junction in an 
environment with sufficiently high impedance is written 
as 

H=^-Ejco S cj>, (1) 

where Q is the charge on the junction electrode, C is the 
capacitance of the junction, Ej is the Josephson energy, 
and <f> is the Josephson-phase difference across the junc- 
tion. The charge Q and hcf>/2e arc quantum-mechanically 
conjugate valuables, and a set of the eigenfunctions are 
Bloch waves of the form 
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FIG. 1: (a) Energy diagram (energy eigenvalue in units of 
Ec vs. quasicharge) and (b) theoretical current-voltage char- 
acteristics for a single Josephson junction with Ej /Ec = 0.3 
and R qp = 10 2 (h/iv 2 e 2 ) at k B T/E c = 0.3, where Ej is the 
Josephson energy, Ec = e 2 /2C is the charging energy, 7? qp is 
the quasiparticle resistance, and ksT is the thermal energy. 
(The energy diagram depends only on Ej/Ec-) 



i>{4>) = u(4>) exp(i0g/2e) , 



(2) 



where q is called quasicharge and u(<p) is a periodic func- 
tion, 



u((j> + 2ir) = u{4>) . 



(3) 



The energy eigenvalue E plotted as a function of q has 
a band structure, and in all the allowed bands, it is 
2e periodic. An example of the energy diagram for 
Ej/Ec = 0.3, where E c = e 2 /2C is the charging energy, 
is shown in Fig. Under constant current bias / x , in 
the absence of quasiparticle or Cooper-pair tunneling, q 
increases uniformly in time according to 



dq 
~dt 



(4) 



so that the state of the system advances toward higher q 
within a given band as time goes on. The average voltage 
is given by 



<U)=^P(z b , (Z ) 



dE(i b ,q) 
dq 



(5) 



where ib is the band index and P(ib, q) is the probability 
that the system is in the state (h,q)- The probability 
P(ib, q) can be calculated by solving a set of coupled 
differential equations of the form 



dP(i b ,q) 
dt 



J2Mib,q,i' b ,q')P(i' b ,q')=0 , (6) 



where the matrix element A(i b , q, i' b , q') describes the rate 
of transition between the states (ib,q) and (i' b ,q'). The 
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dominant process for A(ib, q, i' b , q') depends on the mag- 
nitude of I x . 

An example of the theoretical I-V curve is shown in 
Fig. For sufficiently small I x (region CB), the domi- 
nant process is stochastic quasiparticlc tunneling, where 
q changes by e. This tunneling always occurs along the 
energy parabola q 2 /2C, such that % changes by or +1 
if the initial state is in the lowest band (% = 1) and 
by ±1 for all the other initial states^ The rate for the 
quasiparticle tunneling is given by 



as T — > 0, where 



T(AE) 



AE/e 2 R q 



exp(A£yfc B T) - 1 



(7) 



where R qp is the quasiparticle resistance and AE is the 
difference in energy between the initial (ib,q) and final 
(i' b , q') states, 



AE = E(i' b , q') - E(i b , q) 



(8) 



At sufficiently low temperatures, the tunneling with 
AE > is extremely unfavorable, and the I-V curve 
is highly resistive (Coulomb blockade, region CB in 
Fig. \&). For larger I x (region BO in Fig. [TJ}), the qua- 
sicharge is frequently driven to the boundary of the Bril- 
louin zone, q = e, then taken to — e as a Cooper pair 
tunnels (Bloch oscillation). This process decreases (V), 
and as a result the I-V curve has a region of negative dif- 
ferential resistance, or "back bending," in the low-current 
part. For still larger I x (region ZT in Fig.^Jj), Zcncr tun- 
neling becomes important, and (V) increases again. In 
Zener tunneling, no quasicharge is transferred but the 
state of the system jumps from one band to another as 
it passes by the gap between the bands. The probability 
of Zener tunneling from band i b to % + 1 or vice versa is 
given by 



Pz = exp 



7T (AE) 2 e 



8 i b E c W x 



(9) 



Following Ref.[ll| which takes into account the above tun- 
neling processes (quasiparticle, Cooper pair, and Zener), 
we have calculated the I-V curve numerically. 37 The pa- 
rameters for the calculation are Ej/Ec, ksT/Ec, and 



(10) 



ir 2 e 2 R, 



qp 



The current and the voltage are in units of e/i? qp C and 
e/C, respectively. 

As we have seen in Fig.^i, a typical I-V curve consists 
of three regions, so that it is characterized by the local 
voltage maximum, or blockade voltage V>, and the local 
current minimum, or crossover current I CI . (Here, we 
have to mention that the back-bending feature is blurred 
out if ksT/Ec or a is increased considerably.) An ana- 
lytical expression of V b and I CI has been obtained theo- 
retically for limiting cases^ The value of Vb is a function 
of Ej/Ec, and given by 



V b 



0.25 e/C forSj/^c-Cl 
So/e for EjjEc > 1, 



(11) 



r e 2 / 1 \ 1/4 (Ej\ 3/i 

5o = ^ s {-^) [e^J exp 



Ej_ 

E c 



1/2 



(12) 

is the half- width of the lowest energy-band. As for I CI , 

1/2 



Iz 



RqpC 



(13) 



is expected for a <C (Ej/Ec) 2 <C I and T — ► 0, where 

7T eE? 



Iz = 



hE c 



(14) 



is the Zener breakdown current. Note that I cr is much 
smaller than Iz- When we compare our experimental 
results with the theory, we need theoretical prediction 
for finite ksT/Ec, and arbitrary Ej/Ec and a. For 
this reason we have done the numerical calculation. The 
measured Vb and I CI are compared with the calculation 
in Sec. ESI 



III. EXPERIMENT 
A. Sample preparation and characterization 

A schematic diagram of the AI/AI2O3/AI samples is 
shown in Fig. Et- A single Josephson junction is at the 
center of Fig. |2Jl The junction area is 0.1x0.1 ^m 2 . For 
a scanning electron micrograph of the single junction, 
see Fig. 1 of Ref. 6t On each side of the single junction, 
there are two leads enabling four-point measurements of 
the single junction. A part of each lead close to the single 
junction consists of a ID array of dc SQUIDs. We show 
in Fig- Et> a scanning electron micrograph of a part of the 
SQUID array. Each of the electrodes in the array is con- 
nected to its neighbors by two junctions in parallel, thus 
forming a dc SQUID between nearest neighbors. The 
area of each junction in the SQUID array is 0.3x0.1 /*m 2 
and the effective area of the SQUID loop is 0.7x0.2 /1m 2 . 
All of the samples have the same configuration, except 
for the number N of junction pairs in each SQUID ar- 
ray. The samples are characterized by the normal-state 
resistance R n of the single junction, the normal-state re- 
sistance r' n per junction pair of the SQUID arrays, and 
N. The parameters of the samples are listed in Table [I] 

It is difficult to determine the uniformity of these ID 
SQUID arrays when only the series resistance of all junc- 
tions can be measured. Previous work& showed that the 
variations of the normal-state resistance per junction pair 
were less than 6% for six arrays having nominally iden- 
tical junctions made simultaneously on the same chip, 
with the number of junction pairs ranging from 7 to 255. 
Furthermore, critical-current switching at zero magnetic 
field for these arrays indicated a high degree of uniformity 
(Fig. 4 of Ref. @), and the magnetic-field dependence of 
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TABLE I: List of the samples. R n is the normal-state resistance of the single Josephson junction; Ej is the Josephson energy; 
Ec = e 2 /2C is the charging energy, where the capacitance C is estimated from the junction area assuming 130 fF /^ni 2 ; Vb 
and /„• are the blockade voltage and the crossover current, respectively (see Figs.^i and|SJ; R qp is the quasiparticle resistance, 
which is estimated from the measured J cr in Sec. II V Cl N is the number of junction pairs in each SQUID-array lead; r' n is the 
normal-state resistance per junction pair of the SQUID-array leads. 



Single Josephson junction Leads 



Sample 


Rn (kfi) 


Ej/Ec 


Vb (mV) 


lex (nA) 




N 


r'n (kfi) 


A 


1.0 


9 








255 


1.0 


B 


1.1 


8 








255 


1.1 


C 


1.2 


8 








255 


1.2 


D 


1.8 


5 








65 


0.6 


E 


5.4 


1.7 








33 


0.5 


F 


7.7 


1.2 


8 


0.19 


(0.6 - 1.4) x 10 3 


255 


1.0 


G 


9.0 


1.0 


10 


0.19 


(0.5 - 1.2) x 10 3 


255 


1.1 


H 


12 


0.8 


11 


0.08 


(1 - 3) x 10 3 


33 


1.2 


I 


17 


0.5 


25 


0.08 


(1-3) X 10 3 


65 


1.4 


J 


31 


0.3 


29 


0.02 


(4 - 9) x 10 3 


65 


1.9 


K 


41 


0.21 


28 


0.02 


(4 - 9) x 10 3 


255 


5 


L 


58 


0.15 


29 


0.01 


(6-11) x 10 3 


255 


7 




200 nm 



FIG. 2: (a) Schematic diagram of the samples and the circuit 
for measurements, (b) Scanning electron micrograph of a part 
of the ID SQUID array. 



the arrays was smooth. For the samples of this work, 
the uniformity of the SQUID arrays could be estimated 
by comparing the two arrays on the right- and left- hand 
sides of the single junction. On an average, the normal- 
state resistance of the arrays agreed to within 5%. How- 



ever, if the tunnel-barrier thickness was identical for all 
junctions, one would expect R n /r' n m 6 from the junc- 
tion area for our samples. In Table [IJ we see a large 
variation in R n /r' n {= 1 — 17), which would be mainly 
due to the variation of R n for the following reasons: In 
order to obtain r' n , we do not examine one junction pair 
in the arrays, but measure two arrays in series and di- 
vide by 2N. This kind of "averaging" does not occur 
when we determine R n . Moreover, the single junction is 
much smaller than the junctions in the SQUID arrays, 
and small variations in the lithography and evaporation 
angle cause larger changes in the area of the smaller junc- 
tion, thus leading to larger variations in the junction area 
and thereby the resistance and capacitance. 

The tunnel junctions were fabricated on a SiC^/Si sub- 
strate with electron-beam lithography and a double-angle 
shadow evaporation technique^ We employed a double- 
layer resist (ZEP520/PMGI) method in order to achieve a 
large enough undercut profile for the desired shift, which 
in our case was 0.2 /xm. The electron-gun evaporation 
of Al was done at a rate of 0.1 — 0.4 nm/s in a turbo- 
molecular-pumped vacuum system with the base pres- 
sure of < 2 x 10~ 6 Pa. The vacuum system is equipped 
with an angle stage, and with the stage tilted to an angle 
+9 k 13°, Al was deposited to a thickness of 20 nm. Af- 
ter deposition of the base electrode, O2 was introduced 
into the chamber to a pressure of 1.3 — 4.0 Pa for a period 
of 2 — 10 min. The oxidation conditions, or the thickness 
of the AI2O3 layer, determines R n and r' n . After oxida- 
tion, the top layer of Al was deposited to a thickness of 
30 nm with the stage tilted to an angle —9. 

In the above process, bonding pads (0.3 x 0.35 mm 2 ) 
were also fabricated simultaneously, i.e., the pads are 
20+30=50 nm-thick Al. Onto the pads, 25 /zm-diameter 
Al wires were wedge bonded. 
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B. Measurements 

The samples were measured in a 3 He- 4 He dilution re- 
frigerator (Oxford Instruments, Kelvinox AST-Minisorb) 
at T — 0.02 — 0.6 K. The temperature was determined 
by measuring the resistance of a ruthenium-oxide ther- 
mometer fixed at the mixing chamber. To measure the 
resistance of the thermometer, we used an ac resistance 
bridge (RV-Elekroniikka, AVS-47). Magnetic fields on 
the orders of 1 — 10 mT were applied by means of a su- 
perconducting solenoid. In this magnetic-field range, the 
temperature error of ruthenium-oxide thermometers due 
to the magnetoresistance is negligibly small (e.g., a typi- 
cal value of the error at T — 0.05 K is less than 0.1%). 22 
The samples were placed inside a copper rf-tight box that 
was thermally connected to the mixing chamber. All the 
leads entering the rf-tight box were low-pass filtered by 
1 m of Thermocoax cabled 

We measured the I-V curve of the single junction in a 
four-point configuration (see Fig. EK)- The bias was ap- 
plied through one-pair of SQUID-array leads, and the po- 
tential difference was measured through the other pair of 
SQUID-array leads with a high-input-impedance instru- 
mentation amplifier based on two operational amplifiers 
(Burr-Brown OPA111BM, 10 13 -10 14 fl input impedance 
according to the data sheet) and a preamplifier [Stanford 
Research Systems (SRS) SR560]. The current was mea- 
sured with a current preamplifier (SRS SR570). When 
the voltage drop at the SQUID arrays was much larger 
than that at the single junction, the single junction was 
practically current biased. 

The SQUID arrays could be measured in a two-point 
configuration (same current and voltage leads) on the 
same side of the single junction. Note that the two ar- 
rays are connected in series and that current does not 
flow through the single junction. In order to obtain R n 
and r' n , we measured at T = 2 — 4 K (above the super- 
conducting transition temperature of Al) . 
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FIG. 3: Low-temperature (T = 0.05 K) current-voltage char- 
acteristics of two SQUID-array leads connected in series at 
three magnetic fields for Sample C. The origin of the current 
axis is offset for each curve for clarity. Note the difference in 
the current scale. 



changes depending on B. The I-V curve is Josephson- 
like at B = 0, while a Coulomb blockade is developed at 
B = 7 mT, where the normalized flux 



IV. RESULTS AND DISCUSSION 



/ 



BA 



(16) 



A. SQUID arrays: tunable electromagnetic 
environment for the single junction 

The effective Josephson energy E'j between adjacent 
islands in the SQUID arrays is modulated periodically 
by applying an external magnetic field B perpendicular 
to the SQUID loop, 



E>j<x 



COS 7T 



BA 



(15) 



so long as B is sufficiently smaller than the critical field, 
where A is the effective area of the SQUID loop and 
$o = h/2e = 2 x 10~ 15 Wb is the superconducting flux 
quantum. Figure |3] shows how the I-V curve of two 
SQUID arrays on the same side of the single junction 



is close to 1/2 for our samples with A = 0.14 /xm 2 . Note 
that the current scale for B = is 10 4 times larger than 
that for B = 7 mT. In this manner, the I-V curve is 
varied periodically until B becomes comparable to the 
critical field. At B — 50 mT, the superconductivity is 
suppressed considerably, and the I-V curve is almost lin- 
ear. The magnetic-field dependence of the SQUID arrays 
is summarized in Fig. 01 where the zero-bias resistance R' 
of the same SQUID arrays is plotted vs. B. As expected 
from Eq. 1)150. R' Q oscillates in the low- field regime. In 
the shaded regions, a Coulomb blockade is clearly ob- 
served and the maximum value of R' can be as high as 
10 9 - 10 10 fl (see Fig. 2 of Ref. H for R' of Sample B 
at < B < 9 mT). At B > 40 mT, the oscillation is 
suppressed and R' Q approaches the value in the normal 
state, 2NrL. 
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20 30 
B (mT) 

FIG. 4: Zero-bias resistance of the same SQUID-array leads 
as in Fig. [3] vs. external magnetic field at T = 0.05 K. The 
solid curve is a guide to the eye. The shaded boxes and the 
broken horizontal line represent the region where Coulomb 
blockade was clearly seen in the current-voltage curve and the 
normal-state resistance of the two leads, 2Nr' n , respectively. 
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FIG. 5: Current-voltage characteristics of the single Joseph- 
son junction at T = 0.04 K and / = for Sample G. For the 
definition of /, see Eq. I|16|l The vertical broken lines represent 
the superconducting-gap voltage ±2Ao/e. 



In summary of this subsection, the SQUID arrays, 
which act as an electromagnetic environment for the sin- 
gle junction in this work, can be tuned over a wide range 
in situ with /. Henceforth, we use / defined by Eq. 116(1 
rather than B for the purpose of indicating the external 
field. 



B. Coulomb blockade induced in the single 
Josephson junction by tuning the environment 

In this subsection, we take Sample G for example. 
We look at the I-V curve and the influence of the elec- 
tromagnetic environment on it. (The I-V curves of 
Sample I are shown in Figs. 3-5 of Ref. |y.) Figure [5] 
shows the I-V curves for the single junction at / = 0. 



The vertical broken lines represent the superconducting- 
gap voltage ±2An/e, which is ±0.34 mV for Al ac- 
cording to Ref. 24. In a conventional Josephson junc- 
tion (Ej 3> Ec), a finite current flows at V — as 
expected within the RCSJ (resistively and capacitively 
shunted junction) model. 25 This supercurrent has also 
been observed in small-capacitance Josephson junctions 
[Ej ~ Eq), and when a small-capacitance junction was 
biased with a carefully filtered low-impedance circuit^ 
the magnitude of the supercurrent reached the value pre- 
dicted by the theory. Our single Josephson junction is bi- 
ased with SQUID arrays, whose R' Q is much higher than 
the impedance (~ 24 fi) of the biasing circuit of Ref. |2(] 
even at / = 0. For the case of Fig. [S] (Sample G at 
/ = 0), R' = 0.21 kfl. Ingold and Graberti 5 . consid- 
ered how phase fluctuations suppress the supercurrent in 
a small-capacitance Josephson junction. The small ca- 
pacitance of a junction is a source of phase fluctuations 
because the charge on the capacitance is the conjugate 
variable to the Josephson-phase difference. In addition, 
there is a strong relationship between the fluctuations 
and the electromagnetic environment. They assumed a 
purely Ohmic environment R, and showed that the super- 
current is rounded as R is increased. In such a context, 
the shape of the I-V curve in Fig.[S]may be explained for 
the region in the vicinity of V — 0, although our SQUID 
arrays are not represented by a purely Ohmic impedance. 
At higher voltages, finite currents still flow in Fig.[5]even 
within the superconducting gap. Theoretically, at a fi- 
nite voltage \V\ < 2Ao/e, current can flow when the 
tunneling Cooper pairs can release their excess energy 
2eU, e.g., to the degrees of freedom in the electromag- 
netic environment. Quasiparticle excitations would also 
be responsible for the current especially at \V\ ~ 2Ao/e. 
Below we will concentrate on the immediate vicinity of 
V = (typically ±30 //V), where the I-V curve is most 
sensitive to the state of the electromagnetic environment. 

Figure shows how the single-junction I-V curve de- 
pends on the environment. As / is varied, the curve 
develops a Coulomb blockade. As we have mentioned 
in Sec. [IJ the Josephson energy of the single junction is 
independent of / because it does not have a SQUID con- 
figuration and the field f$>o/A < 7 mT applied here is 
much smaller than the critical field. The electromagnetic 
environment for the single junction (the SQUID arrays), 
however, is strongly varied with /. The behavior of the 
single junction demonstrated in Fig. does not result 
from the magnetic-field influence on the single-junction 
I-V curve, but rather from an environmental effect on 
the single junction. This experiment demonstrates in a 
direct way that the single-junction I-V curve is indeed 
sensitive to the electromagnetic environment. 

Figure0shows the I-V curves of the two SQUID-array 
leads connected in series at the same / as in Fig. |SJ The 
I-V curves of the leads are nonlinear, and in general 
the SQUID array is not described by a liner impedance 
modeL§> However, the environment may be character- 
ized by R' . From (a) to (d) in Fig. |7| R' = 0.21 kfl, 
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FIG. 6: Current-voltage (/-V) curves of the single Josephson 
junction in different environment at T = 0.04 K for Sample G. 
From top to bottom, the normalized flux defined by Eq. 1161 is 
0, 0.40, 0.45, and 0.46, respectively. The origin of the current 
axis is offset for each curve for clarity. 



11 kfi, 0.33 Mfl, and 6.5 MO, respectively. Coulomb 
blockade is distinct only when R' » Rk, which is con- 
sistent with the theoretical conditions for the clear ob- 
servation of Coulomb blockade in single junctions^ For 
an arbitrary linear environment characterized by Z(u>), 
Re[Z(ui)] ^> Rk is required for the Coulomb blockade 
of single-electron tunneling and Re[Z(a;)] » Rq for that 
of Cooper-pair tunnelingiii It is interesting to note that 
while the I-V curve of the leads is "Josephson-like" (dif- 
ferential resistance is lower around V = 0) at all / 
shown in Fig. H that of the single junction is "Coulomb- 
blockade-like" at / = 0.40 and exhibits a clear Coulomb 
blockade at / > 0.45. 



C. Comparison with the numerical calculation 




FIG. 7: Current-voltage curves of the two SQUID-array leads 
connected in series at T — 0.04 K for Sample G. From (a) to 
(d), the normalized flux defined by Eq. HIGH is 0, 0.40, 0.45, 
and 0.46, respectively. Note the difference in the current scale, 
and that the voltage scale for (d) is 10 times larger than that 
for the others. 
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FIG. 8: Blockade voltage Vb divided by e/C as a function of 
Ej/Ec- From top to bottom, the curves represent the numer- 
ical calculations for normalized temperatures ksT/Ec ~ 0, 
0.02, 0.1, and 0.5, respectively. The boxes represent the sam- 
ples with the nominal junction area of 0.01 /jm 2 in Ref. 0- 



The region of negative differential resistance seen at 
the bottom curve of Fig.[fj]is related to coherent tunnel- 
ing of single Cooper pairs according to the theoryi^ of 
a current-biased single Josephson junction in an environ- 
ment with sufficiently high impedance. Following Ref.llll 
we have calculated the blockade voltage VJ, numerically as 
a function of Ej/Ec^ The numerical results are com- 



pared with the measured Vb in Fig. |H1 The data from 
Ref. for the samples with the same nominal junction 
area (0.1 x 0.1 /im 2 ) as in this work are also plotted. 
For Sample G we used the data at / = 0.46 (the bot- 
tom curve in Fig. in order to obtain VJ,. At f — 0.46 
the voltage drop at the SQUID arrays is 10 2 times larger 
than that at the single junction, and the single junction 
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is therefore considered to be current biased. Compare 
the voltage scale of Figs. El and [7|i. We calculated Ej 
from R n , 



For Ec we employed c s = 130 fF //im 2 , and with this 
value the experimental data, especially those for this 
work (the circles in Fig. EJ, agree with the numeri- 
cal calculation. Actually, a smaller value, c s = 45 ± 
5 fF//mi 2 (Ref. l2r|) . which was obtained for the junc- 
tions with 3x28 fj,m 2 and 7 x 54 /im 2 , has been frequently 
employedA^i&il Our apparently large c s may be partly 
explained by distributed capacitance of the SQUID ar- 
rays or by the non-infinite impedance of the environment 
(SQUID arrays), which would have the effect of shunting 
the single junction with an additional capacitance (see 
Fig. Illfl . We also note that the uncertainty in c s seems 
to be large when the junction area is on the order of 
0.01 /jm 2 or smaller, as we have discussed in Ref. |6j. 

Fixing C = 1.3 fF, we tried to reproduce the measured 
I-V curve by optimizing the other parameters of the nu- 
merical calculation, a and k B T / Ec- Examples are shown 
in Fig. EJ where we adjusted a roughly by considering the 
current scale. For both Samples G and I, the slope of the 
experimental data at \V\ < Vb and I ~ can be fit by 
adjusting a, so that it is consistent with the theoretical 
curve for k B T/E c = 0.05 (T = 0.036 K for C = 1.3 fF), 
which is close to the temperature of the mixing chamber, 
0.02 — 0.04 K. At higher bias currents, however, the the- 
ory predicts for ksT/Ec = 0.05 a "back bending" much 
bigger than that in the experimental data. In terms of 
the magnitude of "back bending," the curves for higher 
temperatures, k B T/E c = 0.3 - 0.6 (T = 0.21 - 0.43 K), 
are more similar to the data. "Back bending" in the the- 
oretical curve disappears when the temperature is raised 
to k B T/E c = 0.8 (T = 0.57 K). Although the theory 
explains the general shape of the measured I-V curve, it 
does not yield a complete agreement. It would be rea- 
sonable to assume that the effective temperature of the 
single junction depends on the bias current, and increases 
with increasing current due to heating. This assumption 
reduces the discrepancy between the experimental data 
and the theory. Another possible origin of the discrep- 
ancy would be the imperfect isolation of the single junc- 
tion. While the theory assumes that the single junction 
is completely isolated from the environment, we have bi- 
ased the single junction with the SQUID arrays having 
non-infinite R' . Moreover, the SQUID arrays are a non- 
linear environment. 

Finally in this subsection, we estimate R qp from Id- 
based on the theory. We plot the measured I CI as a 
function of Ej / Ec in Fig. together with some the- 
oretical curves based on our numerical calculation. For 
the estimate of i? qp we used the numerical results for 
k B T/E c = 0.3 and 0.5 (T = 0.21 and 0.36 K), because 
with these values the numerical calculation gave an I-V 
curve similar to the experimental one around I CI as we 




F(uV) 

FIG. 9: Current-voltage curves of single Josephson junc- 
tions, (a) The measured curve (solid circles) of Sample G 
at T = 0.04 K and / = 0.46 is compared with the nu- 
merical calculations for Ej/Ec = 0.53 and a = 0.006 [see 
Eqs. JT0J and JT^J for the definit ion of a and /, respectively], 
(b) Sample I at T = 0.02 K and / = 0.49 with the calcu- 
latins for Ej/Ec = 0.99 and a = 0.003. For both the figures, 
C = 1.3 fF is assumed in the calculations, and from (i) to 
(vi), k B T/E c = 0.05, 0.3, 0.4, 0.5, 0.6, and 0.8, respectively. 



have seen in Fig. El The estimated R qp is listed in Ta- 
ble |U The order of R qp (R qp = 10° - 10 1 MO) and the 
ratio to R n (R qp /R n — 10 2 — 10 3 ) are within a reasonable 
range. 



D. Superconductor-insulator transition in single 
Josephson junctions driven by the environment 

The numerical calculation in the preceding subsection 
agrees with the experiment qualitatively, however, the 
apparently large c s and the discrepancies in Fig. El sug- 
gest that the single junction would still be dependent 
on the environment even at R' ^> Rk- The influence 
of a linear electromagnetic environment Z(u>) on small- 
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FIG. 10: Crossover current I CT vs. Ej/Ec- The curves rep- 
resent the numerical calculations for C = 1.3 fF, and from top 
to bottom (a, k B T/E c ) = (10~ 2 ,0), (ICT 2 , 0.3), (10~ 2 ,0.5), 
(10 _3 ,0), (10~ 3 ,0.3), (10~ 3 ,0.5), (10" 4 ,0), (10" 4 ,0.3), and 
(10 _4 ,0.5), respectively. See Eq. ill IH for the definition of a. 
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FIG. 11: (a) Small-capacitance single Josephson junction 
biased by an ideal voltage source with an electromagnetic en- 
vironment Z(oj) in series, (b) An equivalent circuit to (a) 
when I — V/Z(ll>), where the single junction is biased by an 
ideal current source and shunted by Z(uj). 



capacitance Josephson junction has been studied in the 
framework of the perturbation theoryiilii^ by modeling 
the whole circuit as in Fig. II lb . By assuming an Ohmic 
environment Z(lo) = R, the Cooper-pair current-voltage 
characteristics have been calculated for Ej/Eq <C 1. The 
theory predicts Coulomb blockade of Cooper-pair tunnel- 
ing for Rq / R <C 1 , and supercurrent at or in the vicin- 
ity V = for Rq/R ^> 1. By comparing Fig. |2t with 
Fig. Illt t. we realize that the total impedance Z t seen by 
the junction in Fig. [2^ corresponds to Z in Fig. II lb . If 
we assume in Fig. that the SQUID arrays alone deter- 
mine the effective environment for the single Josephson 
junction at relevant frequencies, where all the four leads 
are shunted at the far ends of the SQUID arrays, and that 
the SQUID array is described by a linear impedance, we 



obtain Z t = (Zf 1 + Z 2 



■ Z x (all 

the SQUID arrays are nominally the same). Note that 
the dc I-V curve is the result of high-frequency qua- 
sicharge dynamics of the single junction, and thus one 
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FIG. 12: Zero-bias resistance Ro of the single Josephson 
junction vs. zero-bias resistance R' of the SQUID-array leads 
at T = 0.04, 0.3, and 0.6 K for Sample G. 



must consider what the junction sees at high frequencies. 
The large capacitance introduced by the cryostat leads 
will effectively shunt the sample with a low impedance 
at the relevant frequencies. This is why we fabricated 
the SQUID arrays in the immediate vicinity of the single 
junction in order to avoid the lead capacitance shunting 
the single junction, and this is why we assume that the 
SQUID arrays alone determine the environment for the 
single junction. We have been characterizing the environ- 
ment by the zero-bias resistance of two arrays measured 
in series, R' . Since Z t ~ Z\, it would be more appro- 
priate to use R' /2. The factor of 2, however, is not very 
important in this work, and thus we keep using R' . We 
should note here that Z ~ R' is a bold approximation. 
The I-V curves of the SQUID arrays are nonlinear as we 
have seen in Figs. |3| and [7| and in general the SQUID ar- 
ray is not described by a liner impedance. The frequency 
dependence 8 is also neglected in the approximation. 



In order to examine the Rq dependence further, we 
plot in Fig. ^| the zero-bias resistance Rq of the single 
junction vs. R' at three different temperatures for Sam- 
ple G. At all the temperatures Ro is a smooth function of 
Rq and increases with increasing Rq, which supports the 
notion that the approximation Z ~ R' is adequate. Even 
at the largest Rq, Rq is still increasing with no sign of 
saturation. The top left point in Fig.[T21for T = 0.04 K 
and Rq = 6.5 Mf2 (/ = 0.46) corresponds to the bot- 
tom I-V curve in Fig. Thus, Rq vs. R' also suggests 
that the single junction would still be influenced by the 
environment at R' 3> Rk- The change of Rq is weak 
at R'q < 10 3 - 10 4 CI. This fact would be understood 
in the following context. In the limit of -Rq — > 0, the 
impedance of the electromagnetic environment for the 
Josephson junction at relevant frequencies is determined 
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by the microwave impedance of the leads, which is usu- 
ally of the order of the free-space impedance Z 377 CI. 
Thus, when R' is of the order of Zq or smaller, our as- 
sumption that the SQUID arrays determine the environ- 
ment is not valid, and Z ~ Zq irrespective of the value of 
Rq. In Ref. 113, special care was taken in order to achieve 
Z/Z < 0.1 at all relevant frequencies. 

The interaction with a dissipative environment, or the 
dissipative dynamics of single Josephson junctions has 
also been studied theoreticallySiiSiil in the circuit shown 
in Fig. Illb with Z(u>) = R. Note that the two circuits 
in Fig. ^] are equivalent when I = V/Z(uj). As the dis- 
sipation is increased, i.e., R is decreased, the Josephson 
junction undergoes a transition from an insulator to a 
superconductor. The phase diagram for this SI transi- 
tion has been derived^ and qualitatively supported by the 
experimentsiiSiiS The phase ("superconductor" or "insu- 
lator") is usually determined by the I-V curve (or dV/dl 
vs. I curved) at the lowest temperatures and/or the tem- 
perature dependence of the zero-bias resistance. 18 We 
also find evidence of the SI transition in the nonlinearity 
of the I-V curve. For Sample G, shown in Fig. H3 we 
see that the I-V curve is "Josephson-like" when the leads 
have R!q = 0.21 kfl (/ = 0), whereas at R' > 0.33 Mfi 
(/ > 0.45), there is a distinct Coulomb blockade. The 
temperature dependence of Rq also supports the exis- 
tence of the SI transition. We see in Fig. ^| that at 
R' < 10 3 — 10 4 ft, Rq has little temperature dependence, 
which is similar to the behavior of a ID systemic in 



the superconducting phase. At R' > 10 4 — 10 5 O, on 
the other hand, we see insulating behavior, where Rq in- 
creases rapidly as the temperature is lowered. 

We show in Fig. ^| the phase diagram for the SI tran- 
sition determined by the I-V curves at the lowest tem- 
peratures of all the samples studied in this work. Note 
that R' was varied over as large as nine orders of magni- 
tude in this work. In Samples F-H, we observed the SI 
transition. The squares and the solid lines to the left of 
the squares indicate that insulating behavior (Coulomb 
blockade or "Coulomb-blockade-like" I-V curve) was ob- 
served. The circles and the solid lines to the right of the 
circles represent the superconducting phase ( "Josephson- 
like" I-V curve). In the range denoted by broken lines, 
the I-V curve was almost linear, and it was hard to de- 
termine the phase. For Samples I-L (Sample E), the 
insulating (superconducting) phase alone was identified. 
Samples A-D have a large Ej/Ec (> 5), and we expect 
from Fig. that their Vb is too small (<C 1 fxV) to be de- 
tected in our dc measurements, i.e., the insulating phase 
is never identified. Indeed, the I-V curve of Samples A- 
D is almost independent of R' . As an example, we show 
in Fig. [H the I-V curve of Sample D at i?' = 0.14 kft 
(/ = 0) and R' Q = 2 Mfl (/ = 0.5). It is Josephson- 
like and in a contrast with the curve in Fig. [3J We do 
not discuss the SI transition for Samples A-D, and just 
indicate the R' Q range covered in the measurements by 
dot-dashed lines in Fig. ^| According to the theory^ of 
single Josephson junctions with an Ohmic shunt R, the 




FIG. 13: Phase diagram for the superconductor-insulator 
transition in single Josephson junctions biased with tunable 
SQUID arrays. The transition is driven by the electromag- 
netic environment characterized by the zero-bias resistance 
Rq of the SQUID arrays. The horizontal axis is the ratio of 
the quantum resistance Rq = h/(2e) 2 ~ 6.5 kfi to R' . The 
vertical axis is the ratio of the Josephson energy Ej to the 
charging energy Ec in the single junction. See text for the 
meaning of the lines and symbols. 
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FIG. 14: Current-voltage characteristics of the single Joseph- 
son junction in two different environment (/ = and 0.5) at 
T = 0.04 K for Sample D. For the definition of /, see Eq. itTol . 



phase boundary is located at Rq/R = 1 irrespective of 
the ratio Ej/Ec when the quasiparticle tunneling is not 
important (a« 1) and the transition is driven solely by 
R. In Samples F-L, a = 10~ 4 - 10~ 2 -C 1 from Table [Q 
It is interesting to see that Fig. ^2 the phase diagram 
for our system with non- Ohmic environment, is consis- 
tent with the theoretical phase diagram where an Ohmic 
shunt is assumed. 
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V. CONCLUSION 

We have studied the current-voltage (/-V) character- 
istics of single Josephson junctions biased with a tun- 
able electromagnetic environment composed of SQUID 
arrays. We have demonstrated that the single junction 
is indeed sensitive to the state of the environment, and a 
Coulomb blockade is induced in the single junction by in- 
creasing the zero-bias resistance R' of the SQUID arrays. 
When R'q is much higher than the quantum resistance 
Rk = h/e 2 k, 26 kfl, a region of negative differential 
resistance has also been observed in the I-V curve of 
the single junction. The negative differential resistance 
is evidence of coherent single- Cooper-pair tunneling ac- 
cording to the theory of current-biased single Josephson 
junctions. Based on the theory we have calculated the I- 
V curves numerically. The calculation has reproduced 
the measured I-V curves at R' 3> Rk qualitatively. 
We have also discussed the superconductor-insulator (SI) 



transition in single Josephson junction driven by the en- 
vironment. We have characterized the environment by 
Rq and determined the phase diagram for the transition. 
This phase diagram is consistent with the theoretical one 
for the SI transition in single Josephson junctions driven 
by an Ohmic shunt alone. 
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